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Abstract 

We study the quantum description of parametric heating in harmonic po- 
tentials both nonperturbatively and perturbatively, having in mind atom traps. 
The first approach establishes an explicit connection between classical and quan- 
tum descriptions; it also gives analytic expressions for properties such as the 
width of fractional frequency parametric resonances. The second approach gives 
an alternative insight into the problem and can be directly extended to take into 
account nonlinear effects. This is specially important for shallow traps. 
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1 Introduction 



Cooling techniques have allowed the trapping and manipulation of atoms by op- 
tical means. Such systems are used to perform experimental tests of fundamen- 
tal principles and have important applications, such as very precise frequency 
standards or studies of cold atomic collisions or collective effects. Fluctuations 
of the electromagnetic fields used to trap or to modify the quantum state lead 
to decoherence in ion traps 0] and limit the trap stability in, e. <?., far-of- 
resonance optical traps As a consequence, there has been, in recent years, 
an increasing interest in understanding the dynamics of noise-induced heating 
in atom traps. 

In the harmonic model of the potential, fluctuations manifest themselves as 
either variations on the spring constant or on the equilibrium position. First 
order perturbative studies [|[ |J| of harmonic parametric heating yield similar 
results to those classically expected: position noise is resonant at the vibrational 
frequency c^o leading to a constant heating rate, while intensity fluctuations 
are resonant at twice the vibrational frequency 2uo leading to an exponential 
growth of the energy. Far-off-resonance optical traps are so sensitive to these 
fluctuations that parametric excitation has been used to accurately measure the 
trap parameters j|, g, ||. 

There are other interesting phenomena predicted by the classical theory of 
harmonic parametric excitation such as resonance effects at fractional frequen- 
cies |Q . They arise from intensity modulations of frequency 2wo / n with n being 
any integer number. The n = 2 resonance has actually been observed in FORT's 
[I, |], Hi an d resonances with n < 10 have been classically studied in ion traps 
§• 

The purpose of this paper is to analyze the quantum description of paramet- 
ric heating in harmonic potentials both nonperturbatively and perturbatively. 
The first treatment shows an explicit connection with the classical problem 
which is valid only for harmonic potentials. It is based on well known algebraic 
techniques |j).It explains observed features of parametric heating predicted both 
classically and quantum mechanically. The second approach has the advantage 
that it can be directly extended to anharmonic potentials. This is specially 
important for shallow traps. 



Nonperturbative treatment of harmonic para- 
metric heating 
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2.1 Time dependent equilibrium position 

The effective Hamiltonian that describes a harmonic oscillator with fluctuations 
in the trap equilibrium position is 



where e q (t) measures such a fluctuation. When the standard creation and 
annihilation a operators are introduced, so that 

9 = \hr— (a + a 1 ) 

V = -i^^(a-at) , (2) 
this Hamiltonian can be written in the form 

H = huJo(a?a + 1/2) + Tif q {t)(a + a 1 ) + %g q (t) (3) 

with 



/«(*) = V ^(*) ! ff«(*) = ^(*)- (4) 
The evolution operator J7(t) satisfies the equation 

dU 

i^=HU (5) 

with the boundary condition 

U(0) = 1. (6) 

Due to the fact that {a^a, at, a, 1} form a closed algebra, it is reasonable to 
write ||1 

Using the well known relation 

e XA Be -XA = B+ * [A B] + ^_ [A [A)B]] + ^ [A} [AB]]] + . . . (8) 

it can be directly shown that the A^'s satisfy the equations 

Ai = -iAaAi +/,(*) 

A 2 = i\ 3 \ 2 + f q (t) 

A3 = 

A 4 = A 3 A 2 Ai-iA 2 Ai+. 9(; (t). (9) 
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with the solution 



ft 

A x (t) = e-' luJot / f q (t')e luJot ' dt' 
Jo 



A 2 (i) = e^ ot / / g (t')e-^ ot df = AJ(t) 
Jo 

A 3 (i) = w i 

A 4 (i) = f (u M\i + i\2(t')Xi(t')+g g (t'))dt', (10) 



which guarantee the boundary condition (|j). We are able, therefore, to evaluate 
the time evolution of any physical quantity. For instance, it turns out that the 
energy operator evolves as 

£(t) = huj Q U- 1 (t)(a^a + l/2)U(t) 

= £(0) + ^o[a 2 a-iAia t + A2Ai)]. (11) 

If e q (t) is a fluctuating held, we are interested in averages over a time T short 
compared with the time scale over which the measurable physical quantities 
change, but large compared with the correlation time of the fluctuations. Let 
us consider fluctuations with zero mean value 

(e q (t)) = ±£ dt'e q (t')=0, (12) 



(e q (t)e(t + a)) = l I dt'e q (t')e q (t' + a) = Vq (\a\). (U) 



and with a stationary correlation function 
Then 

(£(t))=£(0) + ^° I dt' I dt"(e q (t')e q (n)e^''n, , 14 , 



When the one-sided spectrum of the position fluctuations in the trap equilibrium 
position, 

2 f°° 

S q (w) = ~ dacos(ua)i] q (\a\), (15) 
7T Jo 

is introduced and taking t » T, it follows: 



(£(*)) =£(0) + |rn Wo 4 S 9 M, (16) 



recovering the asymptotic expression found by Savard et at However, we 
have obtained it as valid in the nonperturbative regime with only the assumption 
that the fluctuating fields satisfy the equations ( p^HT^ ). Notice that Eq. (O) is 
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also valid whenever the state of the system is a Fock state. This possibility has 
been studied in Ref. Q for the vacuum state where a detailed analysis of the 
short time behavior for different expressions of the correlation function A1A2 is 
performed. 

On the other hand, if one considers driven fluctuations with a well specified 
e q (t) instead of noise fields, the exact equation for the energy evolution (ll|) can 
be used. For instance, if 

e q (t) = e cosuj q t, (17) 

Ai and A2 are trivially calculated and the energy of, say, a coherent state \a) 
evolves as 

(a\£ (t)\a) = hujo(\a\ 2 + -) + ~hLOoeo{ia(C,{{u q - wo)f) + C*(K + wo)t)) - 

- ia*(C*((w, - wo)t) + C(K + "o)t))} + 
t 2 

+ -^o£q[|C(K + ^o)t)\ 2 + |C(K - "o)t)\ 2 + 

+ 2cos(uj t) ^4 2 } (18) 

«y)=e^™f0. (19) 

Equation ( ^8[ ) shows the expected resonances at uj — uj q and emphasizes the 
relevance of the parameter e^a for characterizing the expected heating rates. 

2.2 Time dependent frequency 

In this case the effective Hamiltonian is 



where we have defined 



H=^ + \Mul{l-e{t))q\ (20) 
2m 2 



which takes the form 



H = fo t , o ( t a + 1/2) - ^^l( a + at) 2 . (21) 

Analogously to the former case, the closed algebra nature of {a 1f a + 1/2, a^ 2 , a 2 } 
guarantees that the evolution operator can be written in the form 

U(t) = e co(ata+l/2)/2 e c_atV2 e c + a 2 /2 5 ^ 

whenever the differential equations 

ic -2c_c + = 2u> (l - l/2e(t)) 
i(c_ - 2c 2 _c + )e c ° = -cj e(t) 

ic+e- CQ = \w e(t). (23) 
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arc satisfied together with the boundary conditions 

c (0) = c_(0) =c+(0) =0. (24) 

A connection between the functions {co,c+,c_} and a pair of solutions h±(t) 
and /i 2 (t) of the classical equations of motion 

h k (t)+Lu 2 (l-e(t))h k (t) = (25) 

with the boundary conditions 

hi(0) = l , fei(0) = 

^ 2 (0)=0 , ft2(0) = -£Jo. (26) 



is explicitly given by [O] 



c = -2/n(Mi) 
c+ = M%/Mi 

c_ = -M 2 Mx (27) 



where 



Afi = J[(/n-*Aa) - — (fca + *fti)] 
2 w 

M 2 = h-(h! - i/i 2 ) - — (h 2 + iht)}. (28) 

2 LU 

Eqs. ( gg , g^gg ) show that the theory of classical response of harmonic oscillators 
to noise can be useful in the study of their quantum dynamics. 

Using Eq.(|^), an expression for the time evolution of the energy is found: 

£(t) = huj tf(t)(a'ia+l/2)U(t) 

= hu (l + 2|Af 2 | 2 )[(a t a + 1/2) - MiM 2 a 2 - M£M|a t2 ], (29) 

as well as a nonperturbative expression for nonzero transition probabilities be- 
tween Fock states |fc) and \s): 

um\8)\'= |M2|fc ' s f'f( ir ^ f M) 2 H 2 

I /I 2 fe_ S | Mi | fe+s+ i L^l ; m!(s-2m)!(m+(fc-s)/2)! V 2 7 J 

(30) 

where s and fc have the same parity and for definiteness we have taken k > s. 
Let us focus in the particular case 

e(t) = e cosivt. (31) 
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corresponding to controlled parametric excitation. Then, the classical equations 
of motion can be written in the Mathieu canonical form 

+ (a - 2q cos 2z)h k = 0, (32) 

with z = u>t/2, a — (2ujq/lo) 2 and q — eoa/2. It is well known 12, [ll| that 



depending on the values of a and q, the solutions h k are stable or unstable. 
In the context of parametric heating the transition curves separating regions of 
stability and instability define the width of the corresponding resonance. For 
q << a the resonances are located at a ~ n 2 , i. e., to — 2u>o/n. Their width 
can be found as a power series in eo using, e. g., equations 20.2.25 of Ref. fl2]| . 
Thus, the resonance at 10 = 2loq has a width A2^ ~ eoWo/2 while the resonance 
uj = luq has a width A Wn ~ EqLl>o/6. 

Using the normal form of the unstable solutions of the classical equations of 
motion, it is found that 

Mi = |'e^(t) + ^- ) e-™«/%(-t), (33) 

with 0o a periodic function 0o(i) — <j>o(t + 2tt/u>), and 7 = 7 r + 17^ a complex 
number known as the characteristic exponent. The complex numbers [if are 
determined by the boundary conditions (Eq) and the general form of Mj given 
by Eq. (Bp). Thus, for |7 r w|i >> 1 the energy (|29| ) in the resonance region 
behaves as 

f(t) - ^ [(l + 2|^ +) | 2 elT"l^)(a t a + l/2) + 

- ^ + Vi +) e^>V 2 -^ +) V +) %l^a 2 ]. (34) 

If the initial state is a Fock state, the energy will exhibit an exponential growth 
with a rate determined by the characteristic exponent. 

In order to illustrate these ideas, the time evolution of the transition prob- 
ability between the vacuum and the second excited state is shown in Fig. (1) 
for eo = 0.05, and uj — luq, i. e., the first fractional resonance. We observe 
that for Lot < 1 the transition probability exhibits an approximately polynomial 
growth. However, for longer times an oscillatory behavior arises. In Fig. (2), 
we illustrate the energy evolution starting from vacuum state; we take the same 
values for the parameters to and £q. An exponential growth superposed to the 



oscillatory behavior of function <f)(t), in Eq. (33), is found. The similarities be- 
tween Fig. (lb) and Fig. (2) are due to the fact that the main source of heating 
from vacuum state are precisely the transitions |0) — > |2). These graphs have 
been obtained using a numerical solution of the Mathieu equation. 



3 Perturbative approach to parametric heating. 

The nonperturbative approach to parametric heating made in the last section 
was useful to understand the connection between classical and quantum de- 



7 



scriptions of heating, due to either controlled or stochastic variations of the 
parameters defining a harmonic oscillator. Unfortunately, this approach is valid 
only for quadratic potentials while nonlinear effects may be determinant in ex- 
periments with shallow confining potentials [|| |l4| . The purpose of this section 
is to study some high order perturbative effects due to variations on the strength 
of a confining potential. The resulting equations will be applied to a harmonic 
oscillator but with a straightforward extension to anharmonic potentials. Be- 
sides, this will allow us to understand fractional frequency resonances from an 
alternative point of view. 

The system is described by a Hamiltonian 

H=^- + V(x){l + e v (t)) e v {t)«l. (35) 
Zm 

Following standard time dependent perturbation theory we define the unper- 
turbed Hamiltonian 

H = |- + V(x) (36) 

and work in the interaction picture in which the equation of motion of the state 
is 

ih^^ = e v (t)V(x,t)Mt)). (37) 

The transformed state |^(t)) in the interaction picture is obtained from the 
Schrodinger picture state vector by a time-dependent unitary operator 

\V(t))=e mot / R \*(t)), (38) 

while the interaction operator V(x) is given by 

V(x, t) = e lHotlh V{x)e- lHot/h . (39) 

In this picture the evolution operator satisfies the integral equation 

if 1 

U(t) = 1 - - / e v (t)V(x,t)U(t')dt'. (40) 
" Jo 

Let us consider the transition probability amplitude (k\U(t)\s) between given 
eigenstates of the unperturbed Hamiltonian 

H \n) = E n \n). (41) 

An iterative treatment of Eq. ( p!0| ) gives 

(k\U(t)\s) = 8 ks - l -V ks J^ dt'e v {ty^ *' 

- 72 E V knVns ! dt'evit'y^' [ dfe^t'y^" + ■ (42) 
n „ Jq Jq 
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with ujkn = (Ek — E n )/h and 

V kn =: (k\V\n) = E k 6 kn - (fc||-|n). (43) 
If heating is induced by a controlled modulation of the confining potential 

e v (t) — e cos cut, (44) 

then up to second order in e : 

m + 





1 

= 5 ks + ^eo*Vfc s [C((wfe s - 




+ C((^cs - 


n 


Vnsh — -(C(K S + 

l l(ld n s +0J) 


2w)t) 


+ C,(u) ks t) 




~C((Wfen - 


fu)t) 


- C((^fcn 




■1 x (C((^fe s 

i(w„ s - w) 


2w)t) 


+ ((uJkst) 




e((^fcn - 


+ w)t) 


- C(( w fc« 



)*)) 



These expressions have physical meaning only if the changes in the wave function 
induced by U(t) are small in the interval {0,t). 

For a harmonic oscillator with frequency u>q, one finds that 

hujQ 



V kn = [{2k + l)6 kn + y/k(k - l)5 fe ,„+2 + v / (fc + l)(fc + 2)4,„- 2 ] (46) 

As a consequence, the following products Vk n V ns may be different from zero: 

(i) Vfe ! fe±2Vfc±2,fe±4- In (jig), the resonant terms appear in the combination 
C((4wo ± 2w)t)) — C((2^o i ^t) so that this transition is highly suppressed. 

(ii) VknVnk- Resonances are located at u = 0, 2loq. For to ~ 2wo the contri- 
bution of the transition amplitude is of the form 

2^^M+2| ( ) + \Vk,k-2\ ( )]• 

(47) 

(iii) Vfc,fc±2Vfe±2,fe±2 and Vk,kVh,h±2- These are transitions that may be 
viewed as a combination of two virtual transitions k — > s — > s and fe — > fe — > s. 
The corresponding resonance frequency according to Eq. ([IB]) is the fractional 
frequency a; = |a;fc s |/2 = loq. In fact the transition probability for a modulating 
frequency u> ~ cjo is 

IW (9) Wl a - (^) 4 4i^i a (^-y..) a T a( ""^M 

v 2?i (w — wq) t /4 



^ff(£ - s) 2 [fc(fc - l)£ fc , s+ 2 + (k + l)(k + 2)4 
sin 2 (a; — u )t/2 



8-2 



(48) 
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In all cases, the nonresonant terms ((co't), with u>' ^ 0, give rise to an 
oscillatory behavior of the transition probability which is consistent with that 
found in Fig. 1 for the exact evolution of the transition probability |(0|£/|2)| 2 . 
If one considers sufficiently long times, uj't >> 1, having in mind the delta- 
function representation 

6(w) = - lim ; 1 ' , 49 

it is clear that just the resonant terms have a significant contribution. In such 
a limit the transition probability rates R s <—k are constant. 

Now, some of the general behavior of higher order corrections can be inferred. 
Thus, the dominant transition probability of a fractional frequency resonance 
u> = 2coo/n arises at n-order perturbation theory. It can be interpreted as an 
n-step procedure consisting on n-virtual transitions , where n — 1 of them do 
not change the state and one changes it. Thus, we expect the expression ( fl8[ ) to 
describe approximately the transition probabilities k — > k ± 2 when the source 
has a frequency u> = luq. This can be verified by comparing with the exact 
results of last section. 



4 Discussion and Outlook. 

In this work, we have performed a perturbative and a nonperturbative analysis 
of quantum parametric oscillators. The first approach is based on standard al- 
gebraic techniques and gives a direct connection between classical and quantum 
results. In the case of controlled driving terms of the form e(t) = eocosojt, 
the analytic solutions were used to evaluate time-dependent observables such as 
the energy growth due to parametric heating. This is specially important for 
far-off-resonance-traps (FORT's) when the harmonic oscillator approximation 
is valid. In that case, parametric heating is used as a technique to measure 
the characteristic frequency of the trap. In such experiments, the fractional 
frequency resonances are usually observed and it is clear, within our formal- 
ism, how to perform their quantum description. In particular, using well known 
results of the theory of Mathieu functions, it was shown how to evaluate the frac- 
tional resonances width and how obtain explicit expressions for the exponential 
growth of the energy. 

The possibility of using the nonperturbative analysis to describe noise heat- 
ing effects is a subject that deserves more analysis. This idea has recently been 
exploited JhJ in the case of fluctuations in the trap center. In the case of inten- 
sity fluctuations, the problem is more complicated since a rigorous description 
would require to study the correlation functions of the classical solutions {hi, /12} 
and their derivatives {hi, /12} as they appear in the expressions for the time evo- 
lution of a given physical quantity. For instance, to study the time-evolution 
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of the energy, it would be necessary to evaluate the classical correlation func- 
tions I-M2 1 2 and M\Mi. These correlation functions would depend on the noise 
correlations. 

An alternative would correspond to use solutions {h%, /12} of classical equa- 
tions describing in an effective way the coupling between the harmonic oscillator 
and the fluctuating fields. The usefulness of such an approach would be condi- 
tioned by its ability to reproduce experimental effects. For instance, classical 
harmonic oscillators with an intensity variation e(t) — egCOSLo't and subject to 
a damping force — 7.T lead to fractional frequency resonances only if eo is greater 
than a certain threshold which depends on 7 and the order of the resonance. 
Such thresholds have been observed in the classical collective motion of ions in 
Paul traps ft would be expected that a similar phenomena occurs in the 
quantum regime of motion. The experimental study of these thresholds could 
be used to evaluate effective damping effects in atom traps. 

In this paper, the perturbative analysis of heating induced by variations of 
the intensity of a potential was studied in detail for a harmonic oscillator. This 
analysis gave a different insight into the problem. Fractional frequencies of order 
n appear in n-th order perturbation theory. From the quantum theory point of 
view, they are a direct consequence of: (i)the fact that the harmonic potential 
has diagonal matrix elements different from zero; (ii) the equidistant 

spectrum of the harmonic oscillator. 

The extension of the perturbative analysis to other confining potentials 
is straightforward. In fact, this approach has already been implemented for 
FORT's with shallow potentials @. 

We thank S. Hacyan, P. Nicola and G. Modugno for stimulating discussions. 
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Figure Captions 

Figure 1. 

Nonperturbative time evolution of the transition probability |(0|[/|2)| 2 due 
to parametric heating with a time dependence e(t) = e coswi, when eo = 0.05 
and the first fractional resonance condition ui = luq is satisfied. The short time 
behavior is shown in Fig la. 

Figure 2. 

Nonperturbative evolution of the energy of a time dependent harmonic os- 
cillator of natural frequency uj in units of Tiuj . The system is initially in the 
vacuum state and the intensity time dependence is given by e(t) — eo cos cot, 
eo = 0.05, u> — luq, so that the first fractional resonance condition is satisfied. 
Superposed to the exponential growth, the oscillatory behavior of the function 
4>(t) that appears in the Floquet solution to Mathicu equation is clearly mani- 
fested. 
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